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Abstract — Technological advancement has increasingly 
helped solve engineering problems, but the designer who 
does not seek knowledge of the physical behavior of 
materials and mathematical techniques to solve them runs 
the risk of not being able to interpret the results provided 
by structural software. Computational techniques do not 
provide the exact solution of engineering problems, they 
provide approximations to the intervals where the 
solution of the differential equations governing the 
mechanics of the continuum are inserted. The engineer 
who does not know how to idealize the engineering 
problem in his mind will not be able to model the problem 
in the computer. One of the most intriguing and peculiar 
engineering problems is the plastification phenomenon of 
metallic objects. 

Keywords — Steel Bars, physical behaviour of materials, 
plastification. 

I. INTRODUCTION 

Structural analyzes are made in first-order theory in wich 
the equilibrium of the structure is studied in the 
indeslocated position, and in second-order theory where 
equilibrium is studied in the displaced position. The 
stability of the structure can only be studied in second 
order theory. According to Chen (et. Al., 1996) in the 
study of the resistence and stability of structures it is 
necessary to do: 

• Check the inelastic behavior of the materials; 

• Check the flexibility of the connections; 

• Check second order effects due to geometry 
change. 

The consideration of physical and geometric 
nonlinearities are pre-requisites for the beginning of the 
understanding of the plastification effects of steel 
structures. 


II. ANALYSIS TYPES 



Fig.l: Load-displacement characteristics of general 
analysis types. 

The types of structural analyses are the First and Second 
Order Elastic Analyses and First and Second Order 
Inelastic Analyses, as shown in the figure above. 

In order to understand and delineate the plastification 
effects of steel structures, Chen (1) proposes the 
Advanced Structural Analysis, which consists of a 
technique that seeks to introduce into the behavioral 
models of the elements representative of the structure and 
its materials the closest hypotheses of reality and, 
together, sophisticated numerical and iterative procedures 
to estimate the non-linear behavior of these structures. In 
this way, the advanced analysis encompasses the 
nonlinear, physical and geometric analyzes in the studies 
of the structural systems and their members modeling the 
effects of plastification. 

Three lines of thought for computacional algorithms have 
studied the effects of second order as follows: 

1. Second order inelastic analysis with formation of 
plastic hinge joints and ficticious loads National 
Loads Chen (et al., 1996); 

2. Modified Plastic Hinge Joint Method - allows 
smooth degeneration of stiffness; 

3. Distributed plasticity - sliced discretized bar where 
the inelasticity is modeled with the stress state at 
the center of the slice, Lavall 1996. 
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Fig.2: Plastic-hinge lengths in beams with different 
boundary conditions and load CHEN et al. 1996. 


III. ATTRIBUTES NECESSARY FOR THE 
STUDY OF PLATIFICATION - SECOND 
ORDER EFFECTS. 

The study of plastification in steel elements should be 
done by accounting for second order effects, such as loss 
of physical capacity of steel sheets (inertia, area, thermal 
expansion, etc.), and changes in geometry due to the 
effects of forces. These changes during the action of 
forces must be inserted into the model. The following are 
the items that should be observed in the modeling: 

• Initial imperfections due to the curvature of the bars, 
frames and columns out of alignment, displacement 
of the bars. 

• Residual tensions due to manufacturing and 
assembly processes - Welding action. 

• Connection types: flexible, rigid, semi-rigid. 

• Types of cross sections: symmetrical, non- 
symmetric, open or closed profile. 

• Effects of loading: dynamic, variable, repetitive due 
to construction stages, etc. 

• P-A moment: second-order moments due to axial 
force acting through displacements associated with 
member chord rotation. 

• P-8 moment: second-order moments due to axial 
force acting through displacements associated with 
member curvature bending. 


Plastic phase: (unload) 

(<t - <r y ) > 0 ••• D = D e (4) 

(load) 

D = D e ? (5) 



Fig. 3: Coordinate system and degrees of freedom. 


Corrotational local system: 

q a = (<7i <?2 to} ••• q 1 = l c -l r : (6) 

92 f 

q 3 = a b 

Cartesian degrees of freedom: global system 
Pi = { u a v a e a u b v b 9 b ] (7) 

The structural theory used in this work is based on the 
hypothesis attributed to Bernoulli-Euler, which states 
that: "The cross-sections flat and orthogonal to the axis of 
the bar remain flat, indeformable and orthogonal to the 
axis after deformation." The displacements u c and v c of a 
section point in the corrotational system (x c , y c ) are then 
determined by: 

— y r sena .*• (8) 

v c (x, y) = v c (x) — y r (1 — cosa ) (9) 


IY. FORMULATION 

Engineering strains and stress: 


l c -l r M N 

s = A-l = —— = —; a = a N =— (1) 

Ly Ly il y 


a = a 0) •■•£> = — (2) 


o=> 


Elastic phase: (load e unload); 

(a - (J y ) < 0 D = D e = do Ids (3) 


Stiffiness module or angular 
coefficient or slope of the strain- 
stress curve 


e = (1 + u c )seca — 1 — y r a 


( 10 ) 


Equilibrium Equations: PVW 

Pi=Qaq«i (ii) 


Q a jv r oe a dV r (12) 
(13) 


dP 

dp~ Kt 
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ki,j Qa,i D a ,[: Qfi ,j 4" Qa,i H a ,[l Qfij "I" Qa Qa.ij (14) 

K.P = I V r 8 a D8pdV r H afi = I V r G8 a(i dV r (15) 

Constitutive Stiffness Matrix: Physical and Geometric 
K m = B T DB Physical part (16) 

K g = B t HB + Q a G a Geometric part (17) 

( a = Es -> Elastic Regime 
1 a = Ds Elastic — plastic Regime 


V. COMPUTATIONAL ALGORITHM 

In 1996, Lavall proposed a computational algorithm for 
the study of the plasticity distributed in steel plates 
constituent of beam and frames members. This method 
was known as slicing method based on the formulation of 
item 4. 


6.1 Model 1 - Strength of materials: 

PL 

M m aX — ^ — Mpl — Z X * fy 

E = 20000 kN/cm 2 f y = 25kN/cm 2 


P = 


4 * M, 


pi 


= 2341 kN 


1 ' 


Mp 

L = 500 cm 

Fig. 5: Maximum Plastic concentrated Load Beam 



Figure 4 - Computational tool for evaluating plasticity. 

VI. EXEMPLES 

A static foundation pile test was carried out in the city of 
Rio de Janeiro in 2011 in which the reaction beam was a 
W1000x217 (ASTM A36, fy = 250MPa) with 5 meters in 
length. The test load was 1697kN and so it was done. As 
a research objective to plasticize the bar, the load was 
increased to 2001 kN where from this point the lateral 
buckling mechanism started and the test was stopped. 
This reaction system (W1000x217) is studied in four 
different models: 

1. Model 1 - Strength of materials; 

2. Model 2 - Computational Algorithm of 
Nonlinear Analysis Lavall 1996; 

3. Model 3 - SAP2000 vl9; 

4. Model 4 - Ansys 16. 


6.2 Model 2 - Computational Algorithm of Nonlinear 
Analysis Lavall 1996; - Plastic slices: 



Fig.6: Interface Data Entry Lavall 1996 



Fig. 7: Processing Data , Lavall 1996 
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Fig. 8: Results Data, Lavall 1996 



Fig. 9: Plastified slices zone. 


6.3- Model 3 - SAP2000 vl9: 



Fig. 10: Longitudinal stress crxx web beam. 



6.4 Model 4 - Ansys computer program 16. 




Fig. 14: Longitudinal stress oxx web beam. 


VII. RESULTS ANALYSIS 



Real 

Materials 

Lavall 


Test 

Resistence 

1996 

Maximum 


2341 kN 


2001 kN 

2106,9 kN 

Load 




beginning of 


Result until the ninth 

= 

o 

lateral 


increment of load, from there 

a 

> 

u 


the beam breaks due to 

buckling and 


plastification where the limit 

QJ 

interruption 


of the stress is of 25 kN/cm 2 

& 


and the processing is 

o 



of the test 


interrupted. 



SAP2000 vl9 

Ansys 16.0 

Maximum 

2335 kN 

2341 kN 

Load 

m 

Resolution of the problem up to 

Solving the problem up to the 

# o 

the maximum load, however, 

maximum load and the 

ct 

elements with stress above the 

maximum stresses at the top 

> 


and bottom of the beam are 

Jm 

QJ 

m 

yielding are shown, ex. 235kN / 

shown, however the stresses 

pfi 

o 

cm 2 (2350MPa), which does not 

went beyond the ultimate 


occur in a real test. 

tensile stress 550MPa. 


Table 1 - Comments about the results. 
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